MPW1K Performs Much Better than B3LYP in DFT Calculations on Reactions that Proceed by Proton-Coupled Electron Transfer (PCET)

DFT calculations have been performed with the B3LYP and MPW1K functional on the hydrogen atom abstraction reactions of ethenoxyl with ethenol and of phenoxyl with both phenol and alpha-naphthol. Comparison with the results of G3 calculations shows that B3LYP seriously underestimates the barrier heights for the reaction of ethenoxyl with ethenol by both proton-coupled electron transfer (PCET) and hydrogen atom transfer (HAT) mechanisms. The MPW1K functional also underestimates the barrier heights, but by much less than B3LYP. Similarly, comparison with the results of experiments on the reaction of phenoxyl radical with alpha-naphthol indicates that the barrier height for the preferred PCET mechanism is calculated more accurately by MPW1K than by B3LYP. These findings indicate that the MPW1K functional is much better suited than B3LYP for calculations on hydrogen abstraction reactions by both HAT and PCET mechanisms.     

IR laser-induced dehydrohalogenation of halogenated ethanes

Pulsed IR radiation from a CO₂ TEA laser has been shown to cause dehydrohalogenation of CF₃CH₃, CF₂ClCH₃ and C₂H₅Cl under collisionless conditions. IR emission has been observed and is attributed to HF² and CF₂CH₂² from CF₃CH₃, and to HF² and HCl² from CF₂ClCH₃.     

Are 1,5-disubstituted semibullvalenes that have C2v equilibrium geometries necessarily bishomoaromatic?

Unrestricted density functional theory (UB3LYP), CASSCF, and CASPT2 calculations have been employed to compute the relative energies of the C(s) and C(2v) geometries of several 1,5-disubstituted semibullvalenes. Substitution at these positions with R = F, -CH(2)-, or -O- affords semibullvalenes that are predicted to have C(2v) equilibrium geometries. Calculated singlet-triplet energy splittings and the energies of isodesmic reactions are used to assess the amount of bishomoaromatic character at these geometries. The results of these calculations show that employing strain to destabilize the C(s) geometries of semibullvalenes can lead to a significant decrease in the amount of bishomoaromatic stabilization of the C(2v) geometries, due to reduced through-space interaction between the two allyl groups. However, the C(2v) equilibrium geometries of the 1,5-disubstituted semibullvalenes with R = F and -RR- = -O- do benefit from stabilizing through-bond interactions between the two allyl groups. These interactions involve mixing of the bisallyl HOMO with the low-lying C-F or C-O sigma orbital combinations of the same symmetry. In contrast, for -RR- = -CH(2)-, through-bond interactions destabilize the bisallyl HOMO and are predicted to make the ground state of this semibullvalene a triplet.     

Alkali-metal ion catalysis and inhibition in nucleophilic displacement reactions at carbon,phosphorus and sulfur centres. IX. p-Nitrophenyl diphenyl phosphate

We report on the catalytic effects by alkali-metal ions in the ethanolysis of p-nitrophenyl diphenyl phosphate, in continuation of our studies on alkali-metal ion catalysis and inhibition in nucleophilic displacement reactions at carbon, phosphorus and sulfur centres. The following selectivity order of catalytic reactivity was observed for nucleophilic displacement at the phosphorus center with p-nitrophenoxide as leaving group: Li+ > Na+ > K+ > Cs+. A minor reaction pathway with phenoxide leaving was also found. The kobs data have been dissected into reaction pathways by free ions (kEtO) and by ion pairs (kMOEt), with the latter being dominant, in a 4-membered transition-state. Further analysis is given in terms of initial-state and transition-state stabilization by the alkali-metal ions in terms of the Eisenman model (electrostatic interaction vs. desolvation). Results of ab-initio MO calculations are presented based on interaction between M+ and a model bipyramidal phosphorane intermediate and compared with the sulfurane analogue.     

Ligand-assisted reduction of osmium tetroxide with molecular hydrogen via a [3+2] mechanism

Osmium tetroxide is reduced by molecular hydrogen in the presence of ligands in both polar and nonpolar solvents. In CHCl3 containing pyridine (py) or 1,10-phenanthroline (phen), OsO4 is reduced by H2 to the known Os(VI) dimers L2Os(O)2(mu-O)2Os(O)2L2 (L2 = py2, phen). However, in the absence of ligands in CHCl3 and other nonpolar solvents, OsO4 is unreactive toward H2 over a week at ambient temperatures. In basic aqueous media, H2 reduces OsO4(OH)n(n-) (n = 0, 1, 2) to the isolable Os(VI) complex, OsO2(OH)4(2-), at rates close to that found in py/CHCl3. Depending on the pH, the aqueous reactions are exergonic by deltaG = -20 to -27 kcal mol(-1), based on electrochemical data. The second-order rate constants for the aqueous reactions are larger as the number of coordinated hydroxide ligands increases, k(OsO4) = 1.6(2) x 10(-2) M(-1) s(-1) < k(OsO4(OH)-) = 3.8(4) x 10(-2) M(-1) s(-1) < k(OsO4(OH)2(2-)) = 3.8(4) x 10(-1) M(-1) s(-1). The observation of primary deuterium kinetic isotope effects, k(H2)/k(D2) = 3.1(3) for OsO4 and 3.6(4) for OsO4(OH)-, indicates that the rate-determining step in each case involves H-H bond cleavage. Density functional calculations and thermochemical arguments favor a concerted [3+2] addition of H2 across two oxo groups of OsO4(L)n and argue against H* or H- abstraction from H2 or [2+2] addition of H2 across one Os=O bond. The [3+2] mechanism is analogous to that of alkene addition to OsO4(L)n to form diolates, for which acceleration by added ligands has been extensively documented. The observation that ligands also accelerate H2 addition to OsO4(L)n highlights the analogy between these two reactions.     

Structure and torsional flexibility of the linkage between guanine and fluorene residues in the deoxyguanosine–aminofluorene and deoxyguanosine–acetylaminofluorene carcinogenic adducts

Potential energy surfaces for rotations around two central CN bonds in N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (AAF–dG) and its deacetylated derivative (AF–dG) were studied using Amber 95 molecular mechanics. Both of these adducts are known to be strong mutagens and carcinogens. New Amber 95 force field parameters were derived for the linkage connecting guanine and fluorene moieties in AAF–dG and AF–dG. For this purpose, we determined ab initio MP2/cc-pVDZ//B3-LYP/6-31G* and polarized continuum model Hartree–Fock/6-31G* potential energy surfaces of smaller model systems that included the N-methylimidazole–acetylaniline and N-methylimidazole–aniline adducts. The molecular mechanics parameters were adjusted to minimize differences between the gas-phase ab initio and molecular mechanics surfaces of these model systems. The resulting parameters were transferred to AF–dG and AAF–dG. The barrier for the rotation of the fluorene residue in AF–dG was found to be less than 2 kcal/mol. Such a small barrier renders the fluorene moiety freely rotatable at room temperature. In contrast, the fluorene rotation in AAF–dG is hindered by a significantly larger barrier of 10 kcal/mol. This barrier corresponds to conformations in which the fluorene and acetyl groups lie in the same plane, and is largely due to steric repulsion. Similarly, the coplanar arrangement of guanine and the bridging amino or acetyl groups is disfavored by 5–10 kcal/mol, with AAF–dG again being the more rigid of the two molecules. Energy minima for a rotation around a bond between guanine and the bridging nitrogen are found at ±80° in AAF–dG, and at 120° and –90° for AF–dG. Overall, the fluorene–dG linkages in AF–dG and AAF–dG adducts have significantly different equilibrium structures and torsional flexibilities. These differences may be contributing factors for the observed disparity in mutagenic effects of these adducts.Electronic Supplementary Material: Supplementary material is available in the online version of this article at Acknowledgements. This work was supported by the NSF REU grant no. CHE-0243825 to Loyola University Chicago. We thank to Tom Ellenberger and Shuchismita Dutta for providing us with their results prior to publication.     

Controlling the DNA binding specificity of bHLH proteins through intramolecular interactions

Reversible control of the conformation of proteins was employed to probe the relationship between flexibility and specificity of the basic helix-loop-helix protein MyoD. A fusion protein (apaMyoD) was designed where the basic DNA binding helix of MyoD was stablized by an amino-terminal extension with a sequence derived from the bee venom peptide apamin. The disulfide-stabilized helix from apamin served as a nucleus for a helix that extended for a further ten residues, thereby holding apaMyoD's DNA recognition helix in a predominantly alpha-helical conformation. The thermal stability of the DNA complexes of apaMyoD was increased by 13 degrees C relative to MyoD-bHLH. Measurements of the fluorescence anisotropy change on DNA binding indicated that apaMyoD bound to E-box-containing DNA sequences with enhanced affinity relative to MyoD-bHLH. Consequently, the DNA binding specificity of apaMyoD was increased 10-fold.     

Surface phase behavior and microstructure of lipid/PEG-emulsifier monolayer-coated microbubbles

Langmuir trough methods and fluorescence microscopy were combined to investigate the phase behavior and microstructure of monolayer shells coating micron-scale bubbles (microbubbles) typically used in biomedical applications. The monolayer shell consisted of a homologous series of saturated acyl chain phospholipids and an emulsifier containing a single hydrophobic stearate chain and polyethylene glycol (PEG) head group. PEG-emulsifier was fully miscible with expanded phase lipids and phase separated from condensed phase lipids. Phase coexistence was observed in the form of dark condensed phase lipid domains surrounded by a sea of bright, emulsifier-rich expanded phase. A rich assortment of condensed phase area fractions and domain morphologies, including networks and other novel structures, were observed in each batch of microbubbles. Network domains were reproduced in Langmuir monolayers under conditions of heating–cooling followed by compression–expansion, as well as in microbubble shells that underwent surface flow with slight compression. Domain size decreased with increased cooling rate through the phase transition temperature, and domain branching increased with lipid acyl chain length at high cooling rates. Squeeze-out of the emulsifier at a surface pressure near 35 mN/m was indicated by a plateau in Langmuir isotherms and directly visualized with fluorescence microscopy, although collapse of the solid lipid domains occurred at much higher surface pressures. Compression of the monolayer past the PEG-emulsifier squeeze-out surface pressure resulted in a dark shell composed entirely of lipid. Under certain conditions, the PEG-emulsifier was reincorporated upon subsequent expansion. Factors that affect shell formation and evolution, as well as implications for the rational design of microbubbles in medical applications, are discussed.